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INTERFERENCE HEATING FROM INTERACTIONS OF SHOCKWAVES 
WITH TURBULENT BOUNDARY LAYERS AT MACH 6 * 

By Charles B. Johnson and Louis G. Kaufman n** 

Langley Research Center 

SUMMARY 

This paper presents results of an experimental investigation of interference heat- 
ing resulting from interactions of shock waves and turbulent boundary layers. Pressure 
and heat -transfer distributions were measured on a flat plate in the free stream and on 
the wall of the test section of the Langley Mach 6 high Reynolds number tunnel for 
Reynolds numbers ranging from 2 x 10^ to 400 x 10®, Various incident shock strengths 
were obtained by varying a wedge -shock generator angle (from 10® to 15®) and by placing 
a spherical-shock generator at different vertical positions above the instrumented flat 
plate and tunnel wall. 

The largest heating-rate amplification factors obtained for completely turbulent 
boundary layers were 22,1 for the flat plate and 11.6 for the tunnel wall experiments. 
Maximum heating correlated with peak pressures using a power law with a 0.85 exponent. 
Measured pressure distributions were compared with those calculated using turbulent 
free -interaction pressure rise theories, and separation lengths were compared with 
values calculated by using different methods. 

INTRODUCTION 

In hypersonic flight, the generation of shock waves by various surfaces of a vehicle 
or engine and the impingement of those shocks on other surfaces can greatly amplify the 
local heat transfer and pressure loads on the surfaces. The adverse pressure gradient 
associated with the impinging shock wave, if sufficiently strong, causes the boundary- 
layer flow to separate from the surface. There results a region of reverse flow that 
terminates where the separated boundary layer reattaches to the surface, with attendant 
high heating rates. The importance of being able to predict when boundary-layer sepa- 
ration occurs, the geometric features of the shock— boundary -layer interaction region, 

*This work was Initiated jointly by NASA Marshall Space Flight and Langley 
Research Centers. 

**Staff Scientist, Research Department, Grumman Aerospace Corporation; 
work performed under Langley/Industry Research Associate Program. 
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and the local heating rates is well recognized and has received much attention. (See 
refs. 1 and 2.) However, there are still no adequate analytical methods for predicting the 
extent of separation and the increased heating. (See refs. 3 and 4.) Most of the current 
approaches to the problem of predicting the characteristics of the shock— boundary -layer 
interaction region, particularly for turbulent boundary layers, are empirical and are 
based on data obtained from experiments conducted over relatively limited ranges of 
test conditions. 

Considerable experimental effort has been expended in studying separated flows in 
a shock— boundary -layer Interaction region, a great deal of attention being given to devel- 
oping a method of predicting incipient separation. (See refs* 5 to 12*) Other character- 
istics of separated flows that have been investigated are the extent of the separated flow 
region, the pressure distribution on the surface, and the heat -transfer distribution. 

The extent of separation has been shown experimentally (refs. 13 and 14) to scale 
with boundary -layer thickness and increases with shock strength. Earlier work, based on 
limited range of test conditions (model length Reynolds numbers less than 6 x 10®), indi- 
cated that the extent of turbulent boundary -layer separation increased slightly with 
Reynolds number. (See refs. 13 to 16.) However, Law’s (ref. 6) recent work for larger 
Reynolds numbers indicates a decrease in separation length with increasing Reynolds 
number. The surface pressure rise up to separation can be estimated by assuming a 
’’free interaction.” (See refs. 17 and 18.) The pressure rise downstream of the separa- 
tion point depends on the impinging shock strength and length of separation. (See ref. 19.) 
K the pi-essure distribution is known, then boundary - layer -type analyses can be used to 
estimate the distribution of heating rates. (See ref. 3.) However, as noted by Gerhart 
(ref. 3), the boundary -layer analyses fail at the reattachment point. In the immediate 
vicinity of reattachment, the separated shear layer impinges on the surface, similar to a 
jet impingement type flow, and the heating is most severe. Many simple correlations 
(refs. 15 and 20 to 27) have been sought relating the observed peak heating rates and peak 
pressure values at re attachment. Unfortunately, these correlations are based on data 
obtained for Reynolds numbers considerably smaller (less than 10 x 10®) than those 
expected for high-speed flight vehicles (greater than 50 x 10®). 


Relatively few experiments have been conducted for very high Reynolds numbers 
(refs. 6, 7, 12, and 28), and these investigations have been directed primarily to obtaining 
incipient separation conditions rather than heat -transfer data. The present investigation 
was designed to provide surface pressure and heating-rate distributions resulting from 
incident shock wave impingements for both moderate and very high Reynolds numbers. 
Therefore, an important objective of the investigation was to enlarge substantially the 
base for the empirical methods by providing Interaction <k.ta for a wider range of test 
conditions, and, in particular, to obtain data for Reynolds numbers comparable to those 
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expected for flight vehicles. By using both a flat plate in the free -stream flow and an 
instrumented portion of the tunnel wall, interaction data were obtained for running length 
Reynolds numbers from 2 x 10® to 400 x 10® for unit Reynolds numbers from 9 x 10® to 
120 X 10® per meter. The data were used to improve correlations of separation length 
and peak heating rates and to determine scaling effects on the correlations. The exper- 
iments were conducted in the Langley Mach 6 high Reynolds number tunnel. 

SYMBOLS 

a_ height of sphere -cylinder above plate (tig. 2), m 

Ug height of sphere -cylinder alcove tunnel wall (fig. 2), m 

height of wedge leading edge* above plate (fig. 2), m 
a^ height of wedge leading edge ^ove tunnel wall (fig. 2), m 

Cf Q coefficient of friction of boundary layer on tunnel wall at x = 

* I 

Cm specific heat of model material (stainless steel), 

\ 

Cp specific heat of free-stream flow, JAg-K 

F similarity function used in free -interaction pressure rise (see eq. (8)) 

\ 

f(P) Prandtl-Meyer function of pressure ratio in terms of turning angle 


i 

\ 


^l(Ppk'Ppf) Ppk Ppjg defined in equaUon (10) 

f 2 (ppk>Pp|) function of Ppj. and p^^^ defined in equation (11) 

H ratio of disturbed to undisturbed heat -transfer coefficients at same location 

on plate or tunnel wall 

h heat -transfer coefficient, W/K-m2 

k thermal conductivity, J/sec-m-K 

^sep extent of separated flow upstream of Liviscid, Incident shock location, m 

3 
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Mach number 


Prandtl number 


Stanton number based on undisturbed flow conditions 

exponent in peak heating-pressure correlation power law 
(P/Poo)disturbed 


pressure ratio, 


(P/Poo), 


undisturbed 


pressure, N/m2 

heat -transfer rate, W/m^ 

milt Reynolds number of free -stream flow, per meter 
Reynolds number based on free -stream conditions and distance x 
Reynolds number based on free -stream conditions and distance x 
Reynolds number based on free -stream conditions and distance Xg 
Reynolds mmiber based on free-stream conditions and distance ftg 
recovery factor 

free-interaction pressure rise length coordinate, t 

^ “ ^o 

temperature, K 
time, sec 
velocity, m/sec 

streamwise distance from plate leading edge, m 

distance measured along tunnel wall suiiace from tunnel throat, m 

reference location 
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Subscripts: 


aw 


pk 

pa 


incident shock location on plate, calculated inviscidly, m 

incident shock location on tunnel wall, calculated inviscidly, m 

distance normal to the model surface, m 

ratio of specific heats (taken as 1*4 herein) 

chaise in Prandtl-Meyer turning angle, deg 

boundary -layer thickness, m 

thickness of free-stream layer 

wedge shock generator angle, deg 

viscosity of flow, N-sec/m^ 

density of flow, Icg/m^ 

density of model material (stainless steel), kg/m^ 
thickness of heat -transfer -model wall 
flow deflection angle at reattachment, deg 

adiabatic wall co nditions 

conditions at outer edge of boundary layer 

undisturbed conditions immediately upstream of interaction effects 
peak (maximum) conditions 
plateau conditions 

conditions downstream of reattachment 
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t total (stagnation) conditions 

w wall conditions 

1 undisturbed flow conditions 

free -stream conditions 

APPARATUS AND TEST PROCEDURES 

M 

Test Facility 

The experiments were conducted in the Langley Mach 6 high Reynolds number 
tunnel. This is a closed circular test section, 30.48 cm in diameter, a blowdown wind 
tunnel that can provide Mach 6 free -stream flows for unit Reynolds numbers varying 
from 5 X 10® to 140 x 10^ per meter. (See ref. 29.) 

The tunnel has a model injection system directly beneath the test section (shown in 
fig. 1). Although full injection from the bottom of the vacuum tight chamber to tunnel 
center line requires 1 second, the model passes through the tunnel wall shear layer in 
less than 0.2 second. 


Models 

In order to provide incident shock data for both moderate and high Reynolds num* 
bers, two types of shock receivers were used. Moderate Reynolds numbers were obtained 
on a flat -plate shock receiver mounted in the tunnel free stream. Very high Reynolds 
numbers were achieved by using an instrumented portion of the curved tunnel wall as the 
shock receiver. 

Wedges and spheres were used to generate the shock waves that impinged either 
on the instrumented flat plate or on the instrumented portion of the tunnel wall. For all 
tests, the shock generators v/ere sting mounted to vertical struts. For the “free -stream” 
experiments, the instrumented flat plate was cantilevered forward from the same strut 
(fig. 2). 

The wedge shock generator could be set at 10^, 12.5®, or 15® angles (0), which 
were sufficient to cause moderate extents of separated flow. Smaller wedge angles 
would have resulted in shock strengths insufficient to separate the turbulent boundary- 
layer flows. Larger angles would have resulted in extensive separation and would have 
introduced extraneous effects influencing the flow in the reattachment region, such as 
reflected compression waves and/or the expansion wave from the wedge trailing edge 
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coalescing in the reattachment region. Thus, the attempt was made to provide as large 
a region of "clean” shock impingement interaction flow as possible, without choking the 
tunnel flow. 

As indicated in figure 2, the wedge generator was 19.05 cm long and had a span of 
15.24 cm. The sphere -cylinder shock generator had an overall length of 5.40 cm and a 
diameter of 5.08 cm. Both the wedge and sphere could be set at various streamwise 
stations and at various heights above the shock receivers. The heights are taken in the 
tunnel center plane and are referenced either to the wedge leading edge or to the lower 
surface of the sphere -cylinder. The streamwise stations are referenced to the locations 
where the generator shock waves impinge on the plate surface (xg). The distance Xg 
is measured along the tunnel wall surface from the throat of the tunnel. 

Instrumentation 

Two (interchangeable) flat -plate shock receivers were used, one instrumented 
with 48 thermocouples along the plate center line, and one with pressure orifices at the 
same locations. As indicated in figure 3(a), the plates had chords of 45.60 cm and spans 
of 17.46 cm. The first instrumentation location was 13.94 cm downstream of the leading 
edge, and the thermocouples and pressure orifices v/ere evenly spaced 0.508 cm apart. 
The iron/constantan thermocouples were spotwelded to the inner surface of a "thin skin" 
section of the heat-transfer plate. The skin thickness in this region was 0.076 cm. The 
pressure plate had a constant v/all thickness of 0.318 cm. 

Similarly, two interchangeable tunnel wall sections were used. As indicated in 
figure 3(b), these cylindrical plates were sealed and flush with the tunnel wall. Again, 

48 thermocouples or 48 pressure orifices were evenly spaced (0.508 cm apart) along the 
plate center lines. The first instrumentation location was 3.068 m downstream of the 
tunnel throat (measured along the tunnel wall surface). The heat- transfer plate had a 
thin wall in the region of the thermocouples (0.076 cm thick) , and overall dimensions of 
107.63 cm by 17.46 cm. The radius of curvature of the inner wall was 15.24 cm to con- 
form to the circular test section of the tunnel. The curved pressure plate had the same 
overall dimensions and a wall thickness of 0.318 cm. 

Test Conditions 

The experiments were conducted at five nominal tunnel stagnation pressure levels 
given in table I. Nominal total temperatures and unit Reynolds numbers (per meter) are 
also shown. The corresponding Reynolds numbers based on distance measured along the 
flat -plate surface varied from 1.4 X 10® to 40 x 10®. The Reynolds numbers based on 
distance measured along the tunnel wall from the tunnel throat to the instrumented 
curved plate raried from 20 x 10® to 440 x 10®. 
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Heat-transfer data were obtained for several shock generator locations over the 
full range of test conditions. Pressure data were obtained only for selected shock gen- 
erator locations. 


Test Procedures and Data Reduction 


For both the pressure and heat -transfer experiments, the tunnel flow was first 
started and established, the model injected into the tunnel flow, data recorded, and then 
the model was retracted prior to stopping the turmel flow. Pressures were recorded by 
using standard mercury manometer boards. Heat-transfer data v;ere obtained using 
thermocouples and the thin-wall transient -temperature technique. Schlieren flow photo- 
graphs were taken during each tunnel run, altho .gh in most instances the shock-wave — 
boundary -layer interaction region was outside the tunnel window field of vision. 


The manometers were observed to stabilize within 10 to 15 seconds. However, to 
insure an equilibrium condition, 20 to 30 seconds were allowed before photographing the 
manometer boards. For each pressure run, a second photograph was taken 5 seconds 
later. Readings of both photographs consistently agreed within the accuracy of the 
reading: ±0.1 cm (Ap = ±131 N/ra^ (±0.019 psi)). The resulting uncertainties in terms 
of free-stream static pressures vary from ±0.01p^ for the highest pressure runs to 
±0.19p^ for the lowest pressure runs. 


The temperature time history of each thermocouple was recorded at the rate of 
40 readings/second on magnetic tape by an analog -to -digital data -re cording system. 
These histories were used to obtain the temperature rise rates (dT^/dt), and these were 
used to obtain the heating rates (q^) from 


dT,^ 


( 1 ) 


where and c^ are the density and specific heat of the model material (stainless 
steel), and is the thickness of the thin wall (0.076 cm). Second degree polynomials 
were fitted to the wall temperature readings, taken when the model was positioned in the 
uniform flow of the test section. During injection the model passes through the shear 
layer and enters the tunnel uniform -flow core within 0.2 second. The temperature rise 
rates were taken at the midpoints of the curve fit poiynomjals. 

The error in the heat -transfer data due to longitudinal conduction was checked for 
the case of the spherical generator 1.27 cm above the plate at a unit Reynolds number of 
110 X 10^ per meter. The peak heating for this case has the highest level of heat -transfer 
coefficient of all the data reported and shoxild have the largest conduction error at the 
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point of peak heating. The analysis of ^he conduction first considered the heat balance to 
an element of thin skin which is: 




h(r - T \ 
^ aw *w/ 


Heat stored = Aerodynamic heat in + Heat gained by conduction 

The conduction term on the right-hand side of equation (2) must be evaluated and 
compared with the heat stored term. The second derivative in the conduction term was 
evaluated from the temperature data by a three-point method which would give a maxi- 
mum value of d^T/dx^ at the point of peak heating, and the second derivative was also 
evaluated from second and fourth order polynomial curve fits of the temperature as a 
function of x distance and evaluated at the peak heating. A comparison of the three 
values of d^T/dx^ showed that the three -point method and the second order polynomial 
were approximately the same and gave the largest values. The second derivatives were 
evaluated at various time j and increased with increasing time. The largest value of 
d2T /cbc2 was used to evaluate the conduction term at various times. The results of the 
conduction error check showed that if the heat -transfer data were evaluated at approxi- 
mately 0.3 second after the start of heating, then the maximum conduction error is 
1.7 percent of the indicated heat stored term. For all other peak heating data, the con- 
duction error is less. Therefore, because the conduction error was found to be so small, 
it was neglected in the reduction of the heat -transfer data. 

The experimental heat-transfer coefficient v/as thus calculated from 


T - T 
^aw ^w 


where 


Taw = + r(Tt - T„) (4) 

Because of the prevalence of turbulent boundary layers in these experiments, the turbu- 
lent recovery factor was used in calculating the adiabatic wall temperatuie for the reduc- 
tion of all the heating data; that is, r = Np^ « 0.90 (evaluated at the measured wall 
temperatures). 

The heat-transfer data were obtained during tlie initial heating of the model, while 
the wall was relatively cool ^T^ » 0.6Tt). However, the long time required for the 
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pressure readings to stabilize resulted in the pressures being recorded when the model 
wall was relatively hot « 0.9Tt). Thus, the heating distributions were obtained on 
cool walls, whereas the pressure distributions were obtained on hot walls. Of course, 
this difference in conditions is undesirable, because wall temperature affects both 
boundary-layer transition and separation. This effect must be kept in mind when exam- 
ining and comparing the heating-rate and pressure distributions. 

Interaction Flow Model 

A sketch of the interaction of a shock wave incident to a turbulent boundary layer 
is shown in figure 4, along with typical pressure and heating- rate distributions. In the 
absence of a boundary layer, the incident shock wave compresses and turns the free 
stream flow, strikes, and is reflected from uie wall, and further compresses the stream 
flow. The presence of a boundary layer on the surface complicates tte interaction. The 
incident shock pressure rise causes the boundary layer to separate from the surface 
upstream of the location {kq). where the shock would impinge on the surface if there were 
no boundary layer. As indicated in figure 4, this upstream extent of separated flow is 
referred to as fsep • The pressure rises at separation to a turbulent plateau value, and 
then rises to a peak value in the vicinity where the boundary layer reattaches to the sur- 
face. The boundary -layer thickness is reduced at reattachment. In the vicinity of 
reattachment, there are large peak pressures and very large peak heating rates. 


Flat- Plate - Free Stream 

Undisturbed flow. - Undisturbed (no shock generator) pressure and heating-rate 
distributions on the flat -plate surface were measured for all five nominal tunnel stagna- 
tion pressure levels. These measurements were used as reference conditions to deter- 
mine the pressure rises and heating amplifications caused by the generated shock-v/ave — 
boundary-layer interactions. 

The measured heat -transfer distributions, expressed in terms of Stanton number as 


N 


St,oO 



(5) 


where Cp is the specific heat, is the density, and is the velocity of the free- 
stream flow, are plotted In figure 5 for three Reynolds numbers. The distributions are 
characteristic of transitional -turbulent boundary-layer heating distributions on flat plates, 
and are compared with the theoretical distribution calculated by using an implicit finite - 
difference scheme of Anderson and Lewis (ref. 30) which can be used to calculate the 
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laminai', transitiona?., and turbulent boundary layers. Boundary -layer transition moves 
upstream for the higher total pressure levels, as expected, and at the two highest pres- 
sure levels the boundary layer is characteristically turbulent over approximately two- 
thirds cf the extent of instrumentation on the flat-plate surface (22 cm < x< 39 cm). 

For two of the Reynolds numbers in figure 5, the Anderson -Lewis theory reaches 
the peak turbulent heating at a point slightly upstream of the data, because the theoretical 
location of the end of transition did not match the experimental value. The data and the 
theory both exhibit a peak in heating at the start of fully developed turbulent flow. A 
composite plot of the heat transfer on the flat plate for the four highest Reynolds numbers 
is shown in figure 6 in which a peak at the beginning of turbulent heating is indicated for 
all the da‘‘a. The magnitude of the peak in heating at the start of fully turbulent flow is 
evident from the amount the data are above the Spalding-Chi theory (ref. 31); the calcula- 
tion was made by using a Reynolds number based on distance from the leading edge. 

The difference between the data and theory at peak heating is believed to be a result of the 
theory neglecting the fact that the upstream portion of the flow is made up of a laminar 
and transitional boundary layer. These peaks in heating above flat-plate theory have been 
observed and discussed by Bertram and Neal (ref. 32), and by Neal (ref. 33), and can be 
predicted by a simple flat -plate theory if the correct virtual origin for turbulent flow is 
selected. For a given unit Reynolds number, the data in figure 6 approached the Spalding- 
Chi theory in the downstream part of the fully turbulent boundary layer. 

Shock impingement .- As noted in the previous section, cooling the boimdary layer 
usually delays transition. Thus, transition is believed to have occurred earlier for the 
pressure experiments (T^ « than for the heat-transfer experiments (T^ ^ O.eT^). 

Disturbances, such as the generator shock waves, also cause transition to occur earlier; 
transition in free shear layers (the separated boundary layer) occurs earlier than for 
attached boundary layers (ref. 34). 

In most cases, it is believed the boundary layer was turbulent at reattachment. 
However, for the further upstream shock generator locations and lower tunnel pressure 
levels, the boundary layers were probably transitional in the region of separation. This 
is particularly true for the heat -transfer data, which were obtained during the initial 
heating of the model (cold wall). 

Schlieren flow photographs showing the entire interaction flow region, including 
reattachment, could only be obtained for the upstream shock generator locations, (See 
fig. 7.) The tunnel window fairing obscured the reattachment regions for the downstream 
generator locations. All of the pressure data and much of the heat -transfer data were 
obtained with the shock generators irk their furthest downstream locations. Thus, these 
schlieren photographs do not correspond to the interaction data plotted herein. Indeed, 
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examination of the schlieren photog*'aphs in figure 7 reveals that the extent of separated 
flow decreases with increasing unit Reynolds number for each shock generator. This 
trend indicates transitional separation (boundary -layer transition occurs between separa- 
tion and reattachment). (See refs. 17 and 35.) The thrust herein is toward investigating 
turbulent separation (boundary -layer transition occurring prior to separation); therefore, 
the downstream shock generator locations were used in most cases » It is unfortunate 
that schlieren flow photographs of the interaction region could not be ohi^ed for these 
cases. 

Pressure distributions were obtained only for the furthest downstream locations of 
the shock generators. In these cases, the boundary layer was fully turbulent prior to 
separation at even the smallest free -stream Reynolds number. The measured pressures 
for both the disturbed and undisturbed flat -plate flows were first referenced to the free- 
stream static pressure in order to minimize effects caused by ii^inor variations in p^ 
from run to run, and then the ratios of these pressures were used as the interaction 
pressure distribution: 


P s 




disturbed 


( 6 ) 


(^/^®°)undisturbed 

Wedge shock generators .- Interaction pressure distributions for the 10^ wedge 
shock generator are plotted in figure 8. The unit Reynolds numbers shown are those for 
the tunnel runs with the shock generator installed. (The unit Reynolds numbers for the 
corresponding undisturbed flow tunnel runs were within 10 percent of those shown for tte 
disturbed flow tunnel runs.) The pressure rises abruptly to a value somewhat larger than 
that calculated inviscidly. There is a small additional pressure rise further downstream. 
This rise is attributed to the plate bow wave striking and being reflected by the wedge 
shock generator, and then impinging on the flat plate. This small additional pressure 
rise is followed immediately by the drop in pressure caused by the expansion fan from 
the trailing edge of the shock generator. In figure 8 the inviscid pressure rise and the 
location of shock impingement on the plate surface are indicated by dashed lines. 

The heat -transfer distributions for both the undisturbed flows and for the 10® 
wedge shock interactions are plotted in figure 9. At the lowest Reynolds number 
(R = 17.2 X 10^/m), the undisturbed heating distribution is characteristic of transitional 
boundary- layer flow. However, the heating distributions for the three highest Reynolds 
niunbers are characteristic of fully turbulent boundary-layer flow over the aft half of 
the plate. The interaction heating rate initially drops below the imdisturbed value for 
the lowest Reynolds number case, typical of laminar and transitional separation. No 
similar local regions of decreased heating occur for the three higher Reynolds num- 
ber cases. As in the pressure distributions, there is a small second rise in the 
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heating distributions just upstream of where the heating rates start to decrease. This 
condition is caused by the small additional pressure rise attributed to the reflected 
plate bow wave. 

Distributions of the interaction (disturbed) to undisturbed pressure ratios P, and 
the interaction and undisturbed heat-transfer rate coefficients for the 12.^® and 15^ shock 
generators are plotted in figures 10 to 13. These distributions follow the same general 
trends as described above for the 10^ shock generator. Of course, the pressure and 
heating-rate rises are larger for the stronger shock-wave cases. These distributions 
are for the furthest downstream incident shock location (largest Xg), In every case for 
the three highest tunnel flow unit Reynolds numbers, the boundary layer was fully turbu- 
lent throughout the interaction region. Heating distributions were also obtained for two 
further upstream locations of the wedge shock generators. 

The inviscid, incident shock impingement locations on the flat -plate surface (Xg) 
are shown in table n(a) for the various wedge shock generators. The unit Reynolds num- 
bers (per meter), total temperature, and wall to total temperature ratios are listed for 
each heat-transfer tunnel run. The ratio of the peak heating rate to the undisturbed heat- 
ing rate ^Hp|^ = -pk/^'undisturbed) the corresponding location Xpj^ of the peak heating 
rate are also listed. The pressure tunnel runs were made only for tlie downstream shock 
generator locations. As evidenced in table n, the unit Reynolds numbers for the pressure 
runs approximate those for the corresponding heat -transfer tunnel runs.. However, the 
ratios of wall temperature to total temperature for all pressure runs were approximately 
0.‘J0. The location and value of the peak interaction pressure ratio are listed. For each 
of the three highest Reynolds numbers (fully turbulent separation), the extent (J^sep) of 
the pressure rise upstream of the incident shock location (xg) is indicated. These lengths 
are r- : 2 renced to the calculated (ref. 30) undisturbed boundary -layer thicknesses (6o) at 
the start of the interaction region (these thicknesses are approximately half as large as 
the interval between instrumentation locations on the flat-plate surface). FinaHy, figure 
numbers are listed for the uata plots shown herein. 

Spherical -shock generator .- The interaction pressure ratio distributions for the 
’’higher” spherical-shock generator (ag = 2.54 cm) are plotted in figure 14. In these 
cases there are very abrupt pressure rises to a single measured peak value (for each 
Reynolds number), whereupon the pressure then abruptly falls. Of course, ti^e spherical 
shock wave curves continuously, and the flow passes through a strong expansion fan imme- 
diately downstream of the shock. Because of the very abrupt, spikelike pressure .hstri- 
butions, it is unlikely that the true peak pressure occurred precisely at a pressure tap 
location. Thus, the true peak pressures must be assumed to be somewhat larger than 
those plotted in figure 14 and shown in table H(b) • The inviscid shock impingemem 
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locations Xs were calculated by using an inverse-solution method for sphere-cylinders. 
(See refs. 36 and 37.) 

Plots of the corresponding heating-rate distributions are shown in figure 15. For 
the two lowest Reynolds numbers the boundary layer is transitional throughout the inter- 
action region. The interaction heating distributuions in these cases exhibit the charac- 
teristic initial decrease below the undisturbed values at the start of the interaction. 

There are no local regions of reduced heating for the three highest Reynolds numbers for 
which the boundary layer was fully turbulent throughout the interaction region. 

The strongest shock interaction was generated by the sphere in proximity 
(as = 1.27 cm) to the flat plate. The pressure and heating-rate distributions for these 
interactions (figs. 16 and 17) have the largest gradients and lead to the largest pressure 
and heating-rate amplifications. A second (much smaller) peak is apparent in both the 
pressure and heating-rate distributions. These peaks are associated with a reflected 
shock wave from the sphere -cylinder generator impinging on the flat-plate surface. 


Tunnel Wall 

In order to obtain interaction data for running length Reynolds numbers comparable 
to those expected for flight vehicles, experiments were conducted by using an instru- 
mented section of the wind-tunnel wall as the receiver for the wedge and sphere - 
generated shock waves. In addition to providing thick boundary layers (4 cm < < 6 cm), 

these interaction data provide the extremely wide range of flow conditions necessary to 
verify or correct empirical correlations. 

Undisturbed flow .- Similar to the flat -plate free -stream experiments, undisturbed 
pressure and heating-rate distributions on the tunnel wall were measured for five nominal 
tunnel stagnation pressure levels. These were used as reference conditions for the 
generated shock interaction pressure and heating -rate amplifications. Schlieren flow 
photographs of these Interactions could not be obtained because the tunnel wall is outside 
the window field of vision. 


The variation of boundary -layer thickness on the tunnel wall with Rejmolds number 
is shown in figure 18. The data points are taken from the tunnel wall boundai.7- layer 
surveys of Jones and Feller. (See ref. 29.) The equation faired through these points 
** (see fig. 18) was used to calculate the undisturbed boundary -layer thicknesses used herein. 
Boundary -layer thicknesses calculated by using the theory of Anderson and Lewis (ref. 30), 
with edge Mach numbers and pressures from the method of chai acteristics nozzle design 
program, were approximately 30 percent less than the measured values. The thinner 
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Wedge shock generators As expected, the thick tunnel wall boundary layers sub- 
stantially spread out the pressure rises caused by the interaction of the shock wave with 
the boundary layer, (See fig. 19,) There are ^'plateaus’' in the pressure rises which 
approach the turbulent plateau pressure rise given, for example, by Sterrett and Emery 
(ref. 38) as 


^ = 0.091M? - 0.05 + (7) 

This empirical expression is valid only for Mj > 3.3. The value for Mj - 6, 

- 4.3, is Indicated by the horizontal dashed line in figure 19. 

The incident shock location % is indicated by the vertical dashed line in figure 19. 
The peak pressures occur close to the location of Xg on the tunnel wall. However, 
these peak pressure ratios are substantially less than that calculated inviscidly. (See 
figs. 8 and 19.) The thick tunnel wall boundary layers allow the pressure relief asso- 
ciated with the expansion fan from the wedge trailing edge to propagate forward and 
diminish the peak pressure. The heat-transfer distributions on the tunnel wall, for the 
undisturbed flows as well as for the 10^ wedge shock interaction flows, are plotted in 
figure 20. Oil film flow observations indicated nearly straight separation and reattach- 
ment regions; thus, the transverse curvature effects were negligible. 

The interaction pressure ratio distributions and the heat -transfer coefficient dis- 
tributions for the 12.5^ and 15^ wedge shock generators (figs. 21 to 24) follow the same 
general trends as for the 10® wedge shock generator. Of course, the peak interaction 
pressure ratios, the peak heating amplifications, and the extent of tlie separated flow 
increase with increasing generator shock strength. 

The esqpeirimental results for the wedge-generated shock waves impinging on the 
tunnel wall are summarized in table ni(a). In terms of initial boundary -layer thicknesses, 
the extent of separation ahead of Xg for these tunnel wall tests are comparable to those 
for the free-stream flat -plate experiments for the same peak pressure rises Ppj^. 

(See tables n and m.) For the 15® wedge shock generator, separation occurred upstream 
of the instrumented section of the wind-tunnel wall; this fact precluded obtaining f g^p 
for these cases. 

Spherical shock generator . - Distributions of the interaction pressure ratios and 
heat-transfer coefficients for the spherical -shock wave generator are plotted in figures 25 
to 28. In these cases, the peak pressure ratios, as well as the greatest heating-rate 
amplifications, occurred just upstream of the (inviscidly) calculated shock -Impingement 
locations. For the strongest generated shock wave (% « 2.54 cm, fig. 28), there are 
small secondary pealcs in the interaction heat-transfer coefficient distributions. These 
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peaks are associated with compression waves being reflected from the shock generator 
and interacting with the tunnel wall boundary layer. 

The distance between the tunnel wall and sphere is comparable to the thickness of 
the tunnel wall boundary layer. Therefore, in all these cases, the sphere shock is pre- 
dominantly within tlie rotational flow shear layer of the thick tunnel wall boundary layers. 
This location reduces the shock strength and results in the very small (£gep /^o) values 
shown in table in(b). This condition must be taken into account when comparing the 
experimental results with analytical ones, and when using the experimental results in 
extending or modifying empirical correlations. - ' 

THEORETICAL ANALYSES, COMPARISONS, AND CORRELATIONS 


Free Interaction Pressure Rises 

Although the extent of separation and the total -pressure rise depend strongly on the 
shock strength, the initial pressure rise to the plateau value should be independent of the 
mechanism generating the shock wave. This concept of a ’Tree interaction’* between the 
viscid boundary layer and the inviscid external flow in the upstream part of toe interaction 
region has been well established for a limited range of Reynolds numbers. (See refs. 17, 
18, and 39.) 


Initial pressure rises for the 10^ and 12.5® wedge shocks impinging on the tunnel 
wall boundary layers are indicated in figure 29. The pressure ratios are plotted against 
(x - where the subscript o indicates conditions at the upstream location where 

the disturbance pressure rise begins. The pressure rise curves, thus shifted so that 
they all start at x = Xq, are all very similar. There is no consistent Reynolds number 
effect. However, the 10® wedge shock data indicate a somewhat larger pressure gradient 
than the 12.5® wedge shock data do. Unfortunately, the 15® wedge shock data cannot be 
included in figure 29 because the pressure rises started upstream of the instrumentation 
region (see fig. 23), and therefore the values of Xq could not be determined for these 
cases. 


The effects of the different wedge shock generator angles are eliminated by plotting 
the data in terms of the universal similarity function F(s) suggested by Carriere, 
Sirieix, and Solignac (ref, 18) 


F{s) 


2(P ■ 1) 
yM| 



( 8 ) 
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where s s (x - Xo)/(xf - Xq), Ai/ is the change in the Prandtl-Meyer angle for isentropic 
compression Au = f(P), y is the ratio of specific heats (y = 1.4 herein), c^ ^ is the 
coefficient of friction at Xq, and X| is a reference location. This location, which scales 
the extent of the pressure rise, is determined by requirir^ that F = 4.22 at x = Xf . 

It has been found (ref. 18) that by using this empirical value, turbulent -boundary-layer 
pressure rises correlate well for a wide range of test conditions. 

Jn figure 29, the pressure rises are all forced through P = 1 at x = Xq; in fig- 
ure 30, the pressure rises are all forced through F = 4,22 at x = X|. Although the 
similarity function F does eliminate effects of different wedge shock generator angles, 
there remains a small and inconsistent Reynolds number effect. Included in figure 30 is 
the free interaction pressure rise curve suggested by Erdos and Pallone (ref. 39) for 
turbulent boundary layers. The present data exhibit somewhat steeper initial pressure 
rises than the curves shown in figure 30, although the data for tht intermediate Reynolds 
numbers (Re,x = 211 x 10® and 213 x 10®) fall close to the curves for both wedge angles- 

Elfstrom (ref. 10) and Reeves (ref. 19) observed increases in turbulent plateau 
pressures with increasing separation lengths. This condition holds true for separation 
lengths up to 10 times as long as the undisturbed boundary-layer thickness. The data 
herein support this trend. (See figs, 19, 21, and 23.) However, as noted by Reeves 
(ref. 19) and also by Whitehead, et al. (ref. 4), three-dimensional flow effects are impor-^ 
tant and can alter the trends established for two-dimensional iiows. 


Separation Lengths 

The separation lengths listed in tables II and in, nondimensionalized with respect 
to the undisturbed boundary -layer thicloiesses at the start of the interaction regions, are 
plotted against peak pressure rise ratios in figure 31. These results can be approximated 
very simply by 


e p 

sep _ pk 


( 9 ) 


for Mj = 6 and Reynolds numbers from 10*^ to 4 x lO®. The ^ncn circle symbols in 
figure 31 are for the sphere -generated shock waves impinging on the thick tunnel wall 
boundary layer. In these cases most of the sphere shock is within the rotational flow of 
the boundary layer; this fact makes questionable the ns'j of = 6 as the effective free- 
stream Mach number. With the exception of these points, the remaining data (for these 
particular test conditions) follow equation (9). The data scatter for the sphere -generated 
shock waves impinging on the flat-plate boundary layer can result from the true peak 
pressure not being measured (see figs. 14 and 16) or from lateral relief of the separation 
region. 



I 


I 


Williams (ref. 13) and Mikesell (ref. 14) established expressions for turbulent 
separation lengths in terms of the peak pressure Ppj^ and the turbulent plateau pressure 
and corresponding Mach number Pp£ and Mp£. Williams presented his separation 
length data as a function of the parameter 


whereas Mikesell presented his data as a function of the parameter 


*2(Ppk’Pp«) " 


1 + 


' m2 


nfl 

r 




2 “p4PpJ 


( 10 ) 


( 11 ) 


The separation lengths are functions of these parameters. These functions are shown in 
figure 32 along with the experimental results. Because the plateau conditions must be 
known, only the data for the wedge -generated shock wave impinging on the tunnel wall 
could be shown in this figure. Also, the 15® wedge data had to be excluded because the 
extent of separation is not known for these cases. (See fig. 23.) 


Peak Heating - Pressure Correlations 

The turbulent interference heating correlation shown in figure 33 presents the peak 
heating ratio as a function of the peak pressure ratio. The peak heating ratio Hpj^ is 
defined as the peak heating rate for the interaction flow divided by the heating rate for 
undisturbed flow at the same location. The peak pressure ratio Ppi^ is defined as 
Ppk/Poo disturbed flow divided by p/p^ for the undisturbed flow at the same 

location. The data from the present investigation are shown as various symbols that 
indicate the nominal unit Reynolds number, type of generator, and whether the test is on 
the plate or on the tunnel wall. The curved-line shaded area represents data from vari- 
ous investigators presented by Holden in reference 40. The data appear to correlate 
according to the power law 

Hpk = (Ppk)“ (12) 

The same power law correlation has been presented in many previous investigations. 

‘rhe values of n for the best correlations were found to be n = 0.85 by Holden (ref. 40) 
and Markarian (ref. 26); n = 0.80 by Sayano (ref. 27), Haslett et al. (ref. 22), and Hung 
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and Barnett (ref, 21); and n = 0,71 by Popinsid (ref. 15). The data from the present 
experiments, which include results for ninning length Reynolds numbers exceeding 
400 X 10®, coupled with the previously reported data shown in figure 33, indicate that the 
best correlation is achieved by using n = 0.85. There are no discernible Reynolds 
number effects on the correlation for running length Reynolds numbers varying from 
2 X 10® to 400 X 10®. 


Another correlation, suggested by Bushnell and Weinstein (ref. 25), Is shown in 
figure 34 as 


hpk p. / PwVe®s 

PwVp \'"wsin<^)y 



where PwVe wd /Xvv pertain to the reattached flow, 63 Is the thickness of the free - 
shear layer, <j> Is the turning angle of the flow at reattachment, and subscript e refers 
to conditions at the outer edge of* the boundary layer. The data presented by Bushnell 
and Weinstein were measured at the reattachment of separated flows on trailing -edge 
flaps. Their data exhibit lower heating rates than those obtained from the present exper- 
iments for reattachment of flows separated by incident shock waves. Included in fig- 
ure 34 is a correlation of incident shock wave interaction data presented by Keyes and 
Hains. (See ref. 34.) 


CONCLUDING REMARKS 

Measurements of pressure and heat transfer in the area of a shock-wave— 
boundary -layer interaction were obtained on a flat plate in • free stream and on the 
test-section tunnel wall for Mach 6 and for running- length Reynolds numbers from 2 x 10® 
to 400 X 10®. The test results were obtained for relatively thin (flat plate) as well as 
relatively thick (tunnel wall) boundary layers; the thick boundary layers occur at Reynolds 
numbers comparable with those expected for flight vehicles (running-length Reynolds num- 
bers exceeding 400 x I06). 

The highest interference heating for a turbulent boundary -layer interaction was 
obtained on the flat plate in the presence of a strong shock generated by a sphere -cylinder 
in proximity to the plate. This interaction resulted in amplification factors, compared 
with undisturbed values, of 22.1 and 30.3 for heating rates and pressures, respectively. 
The highest interference heating on the tunnel wall occurred when a shock from a 15° 
wedge generator impinged on a turbulent boundary layer and resulted in disturbed to 
undisturbed amplification factors of 11.6 and 15.4 for heating rates and pressures, 
respectively. For turbulent boundary -layer flow, the peak beating ampLIication can be 
estimated from Hpj^ where Hp^ is the peak ratio of the heat-transfer 


coefficients and Ppj^ is the peak pressure ratio. There are no discernible Reynolds 
number effects for Reynolds numbers from 2 x 10® to 400 x i06. 

The initial pressure rises on the tunnel wall follow the turbulent free interaction 
pressure rise distribution. Free interaction similarity methods adequately account for 
the effect of different shock strengths. However, the data still showed inconsistent 
Reynolds number effects. 

' Separation lengths for these' turbulent boundary layers varied approximately line- 
arly '^th the peak pressure rises. The measured separation lengths are comparable 
with values calculated using previously established empirical expressions, which are 
functions of the turbulent plateau pressure and peak pressure. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., June 3, 1974. 
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TABLE n,- FREE-STREAM FLAT-PLATE EXPERIMENTAL RESULTS 


(a) Wedge shock generator 








1 - 

) 


TABLE a* FREE-STREAM FLAT-PLATE EXPERIMENTAL RESULTS - Concluded 


(b) Spherical shock generator 



*At the lowest tunnel stagnation pressure level, the tunnel flow Mach number approximately equals 5.9, and results In a slightly 
different Inviscld shock shape than for all other tunnel stagnation pressure levels, for which Mj « 6,0. 
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TABLE EDL- TUNNEL Vr ALL EXPERIMENTAL RESULTS - Concluded 


(b) ^herlcal shock generator 



*At the lowest tunnel stagnation pressure level, the tunnel flow Mach number approximately equals 5.9, and results In a slightly different ^viscid 
shock shape than for all other tunnel stagnation pressure levels, for which Mj « 6.0. 
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I^umel wall tests 

Figure 2.- Model sketches and coordinate system. 
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(a) Flat plate. 

Figure 3.- Model instrumentation locations. 
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Figure 6.- Transitional and turbulent heating distributions for various unit Reynolds numbers. 
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(a) 10® wedge generator. 


Figxire 7.- Schlieren flow photographs, generators in forward position, 
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(d) Sphere generator (2.54 cm above plate). 
Figure 7.- Continued. 
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(e) Sphere generator (1.27 cm above plate) 
Figure 7.- Concluded. 
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Figure 9.- Reynolds number effects on the heat -transfer coefficient distributions 
on fist-plate, 10° wedge shock generator. 
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Figure 11.- Reynolds number effects on the heat-transfer coefficient distributions 
on flat-plate, 1?.5^ wedge shock generator. 





Figure 12.- Interaction pressure ratio distributions on flat-plate, 
15° wedge shock generator. 
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Figure 15.- Reynolds number effects on the heat-transfer coefficient distributions 
on flat plate, spherical- shock generator, s 2.54 cm. 
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Figure 16.- Interaction pressure ratio distributions on flat -plate, spherical-shock 

generator, = 1.27 cm. 
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Fl^e 17.- Reynolds ntimber effects on the heat-traixsfer coefficient distributions 
on fiat-plate, spherical- shock generator, a_ * 1.27 cm. 
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Figure 19.- Interaction pressure ratio distributions on tunnel wall 
10® wedge shock generator. 
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Figure 20.- Reynolds number effects on the heat-transfer coefficient distributions 
on the tunnel whll, 10° wedge shock generator. 
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Figure 22.- Reyiiolds number effects on the heat -transfer coefficient distributions 
on the tunnel wall, 12.5® wedge shock generator. 
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Figure 23.- Interaction pressure ratio distributions on tunnel wall, 
15® wedge shock generator. 
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Figure 24.- Reynolds number effects on the heat-transfer coefficient distributions 
on the tunnel wall, 15° wedge shock generator. 


1 

I 

I 


61 




Surface distance from nozzle throat, 5^. m 


Figure 25.- Interaction yressure ratio distributions on tunnel wall, sphericat-shock 

generator, tig » 5.08 cm. 
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Figure 26.- Reynolds number effects on the heat-transfer coefficient distributions 
on the tunnel v/all, spherical-shock generator, Sg * 5.08 cm. 
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Figure 27.- Interaction pressure ratio distributions on tunnel wall, spherical-shock 

generator, Sg = 2.54 cm. 
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Figure 28.- Reynolds number effects on the beat-transfer coefficient distributions 
on the tunnel wall, spherical-shock generator, S- s 2.54 cm. 
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Figure 29.- Initial portion of interaction pressure ratio rises on tunnel wall. 
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Figure 34.- Correlation of peak heating at reattachment for fully turbulent boundary layers, using parameters 

suggested by Bushnell and Weinstein (ref. 25). 
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